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The effective potential of electron-electron interaction and the two-particle “density-density” 
correlation function have been calculated for a simple semiinfinite metal making allowance for 
the local-field correction. The influences of a flat interface and various models of local-field 
correction on the results of calculations are analyzed. 
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1. Introduction 

The modern quantum-mechanical statistical theory 
of Fermi systems with an interface is still far from 
being completed. The urgency in the theoretical de¬ 
scription of such systems can hardly be overestimated 
owing to the importance of processes that occur in the 
presence of an interface and a rapid development of 
experimental methods aimed at surface researches. 

The most popular theoretical method to study such 
systems is the density functional one [H, which was 
created on the basis of the well-known Thomas-Fermi 
approximation firstly developed for atoms. By its na¬ 
ture, the density functional method is a one-particle 
approach, so it cannot consider many-body correla¬ 
tion effects correctly. Therefore, the energy function¬ 
als of systems with an interface are most often consid¬ 
ered in the local density approximation; namely, the 
expressions known from the theory of uniform sys¬ 
tems are taken for calculations, but the distribution 
of the electron density n(r) is substituted for the elec¬ 
tron concentration n. This approach is debatable [3|, 
because the interface introduces not only quantita¬ 
tive, but also qualitative changes in various param- 
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eters of the electron system (e.g., the image forces 
emerge and so on), which cannot be taken into ac¬ 
count in principle by the density functional theory. 

In the cycle of works Ml , an attempt was made to 
develop a consistent quantum-mechanical statistical 
theory for a simple metal with the interface “metal- 
vacuum”. In particular, it was shown that the ther¬ 
modynamic potential and the structural distribution 
functions of electrons in a semiinfinite metal are ex¬ 
pressed in terms of the effective electron-electron in¬ 
teraction potential. This work logically continues this 
cycle of works. It was aimed at studying the influ¬ 
ence of various approximations for the local-field cor¬ 
rection on the two-particle correlation function “den¬ 
sity-density” and the effective potential of electron- 
electron interaction. 

The effective interaction between charged particles 
in spatially confined systems attracts the attention 
of researchers for a l ong time. In particular, it was 
studied in works UMl. In work [i(tI |. under cer¬ 
tain approximations, the polarization part of the en¬ 
ergy of interaction between a motionless point charge 
and a semiinfinite metal was calculated. Using similar 
approximations, the dielectric function of a semiinfi¬ 
nite metal and the effective potential were calculated 
l^ll3| in the framework of 
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in work [ll|. In works 
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the effective interaction potential approximation, the 
asymptotics of Friedel oscillations at large distances 
between charges at the interface was obtained, and 
their dependence on the Fermi surface shape in met¬ 
als were studied. 

In works [l3 - [l^ . the screening of a charged impu¬ 
rity near the metal surface was studied in the random- 
phase approximation. In particular, the dependences 
of the electrostatic potential of this impurity on the 
distance in the surface plane and the distance to the 
surface were calculated in the quasiclassical case and 
the Thomas-Fermi approximation, and Friedel oscil¬ 
lations were revealed. 

In order to eliminate shortcomings inherent to the 
density functional method, attempts were made to 
use the GW approach [3| with the density functional 
method. In particular, in work (l8l |. the effective po¬ 
tential of electron-electron interaction was calculated 
with the use of the local density approximation for 
the exchange-correlation potential. 

2. Model 

Consider a semiiiifiiiite simple metal in the framework 
of the jellium model, i.e. when the ionic subsystem 
of a metal is represented by a uniformly distributed 
positive charge confined by the interface plane z = 
= —d, with the density 

n+{x, y, z) = n+{z) = Ubuik 0{-d - z), (1) 


where 



is the Heaviside function, ribuik is the electron con¬ 
centration, and d > 0 is a parameter determined self- 
consistently from the electroneutrality condition 


+ 00 

J dz (n(z)-n+(z)) = 0, 

— OO 


( 2 ) 


This model of potential physically correctly corre¬ 
sponds to a real situation and allows analytical so¬ 
lutions to be obtained for the Schrodinger equation 


A 

-A ■ 

2m 


V{z) 




£^a(p)^'p.a(r), 


where m is the electron mass, hp the two-dimensional 
vector of electron momentum in the plane parallel 
to the interface, and a the quantum number associ¬ 
ated with the electron motion normally to the inter¬ 
face. Since the electron moves freely in parallel to the 
interface, the wave function and the corresponding 
energy of the electron can be written as follows: 




h^{p^ + a^) 
2m 


Here, r = (r||, z) is the radius vector of the electron, 
S is the interface area, and 


/ \ 2 ./ \/i/ \ 27m 

‘Paiz) = sm(az)0{-z), a = —, n = 1, 2,..., 


where L determines the variation range of the 
electron coordinate normal to the interface: z € 
G [—L/2,-l-oo). The parameters L and S tend to in¬ 
finity, so that the problem is considered in the ther¬ 
modynamic limit. 


3. Effective Potential 
of Electron-Electron Interaction 

According to work Q, the two-dimensional Fourier 
transform of the effective potential of electron- 
electron interaction with respect to the radius vector 
r|| can be expressed in the form 


5(<7ki,Z2) = p(g|zi - Z2) + 

+ 130 +130 


ydzy dzV(q|zi — z) X 


SL^ 

-L/2 -L/2 

X m,{q\z,z')v{q\z' - Z2), 


( 4 ) 


where n{z) is the electron density distribution. Let 
the ionic subsystem form a surface potential for elec¬ 
trons in the metal, which does not allow them to es¬ 
cape. This surface potential is simulated by the po¬ 
tential wall 
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where 

v{q\zi-Z2) = 

q 

is the two-dimensional Fourier transform of the 
Coulomb potential, zi and Z 2 are the coordinates of 
electrons reckoned normally to the interface, /3 the 
inverse thermodynamic temperature, and SJl(q\z,z') 
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the “density-density” correlation function, which is a 
solution of the Fredholm integral equation of the sec¬ 
ond kind 0, 

M{q\zi,Z2) = Tlo{q\zi,Z2) + 


SL^ 


+ 00 +00 


dz Jdz'D}lo(q\zi, z) x 


-L/2 -LI2 


X (^{q\z - z') - iz{q\z - z'jj M{q\z', Z 2 ), 


( 5 ) 


where izkiq) = Gk{q)i^k{q), Vk{q) = 4^ is the 
three-dimensional Fourier transform of the Coulomb 
potential (the variable k is responsible for the Fourier 
expansion along the electron coordinate normal to the 
interface), 


7(g|z-z') = y^e' 


k{z-z’)- 


vk{q), 


d}lt){q\zi, Z 2 ) is the two-particle “density-density” cor¬ 
relation function in the ideal exchange approximation 
0 , and Gfe {q) is the local-field correction. 

In the low-temperature limit, the following expres¬ 
sion was obtained in work @ for the two-particle cor¬ 
relation function of an electron gas in the ideal ex¬ 
change approximation: 

L'i 

MQ{q\zi,Z2) = ^ X 

ai ,Q!2 

X V’l, {Zl)^a2 (^ 2 ), (6) 


where 
^ai,a2 io) 


2m S al — a 2 — q^ 
Tfi 2 tt 



Pf 

(a? - ai - g2)2 
Py - a 


(a? - ai - q^Y 


9{py - ai), 


( 7 ) 


Pf = (97r/4)^/^/rs is the Fermi momentum, and rs 
the Brueckner parameter in the units of the Bohr ra¬ 
dius Ob- 

In work it was shown that, in some approxi¬ 
mations, an analytical expression can be obtained for 
the function ^(}lo(q\zi, Z 2 ), 


mQ{q\z,z') 


SL'^ 2m 1 e-‘?l^-^'l - e-«l^+^'l 
f} h? q 


f pF cos(pf(-z +-z')) _ Pf cos(pf(-z - z')) 
y {z + z'Y {z-z'Y 

sin(pF(z — 2;')) sin(pF( 2 -I-2')) 

+ (z + z'Y . 

xe{-z)9i-z'). (8) 

From this expression, one can see, in particular, 
that, besides the terms, which depend on (z — z') 
and are characteristic of uniform systems, there are 
terms depending on {z + z'). The reason for their ap¬ 
pearance is the presence of the plane interface. This 
means that, in our ideal-exchange approximation for 
the two-particle “density-density” correlation func¬ 
tion, the image-force effects are already taken into 
consideration. Note that the polarization operator in 
works Elm is presented as a sum of two po¬ 
larization operators for the uniform electron gas: 
one of them depends on {z — z'), and the other on 
{z — z'). However, formula ([5]) shows that the depen¬ 
dences of Dyi()(q\zi, Z 2 ) on {z — z') and {z + z') are not 
so simple. 

In the following numerical calculations of the two- 
particle “density-density” correlation function in the 
ideal exchange approximation, expression ([S]) is used. 


4. Results of Numerical Calculations of 
Two-Particle “Density Density” Correlation 
Idmction and Effective Electron-Electron 
Interaction Potential 

The two-particle “density-density” correlation func¬ 
tion of electrons, S[Jl(g|zi, 2 : 2 ), was numerically calcu¬ 
lated according to Eq. ([5]) and making allowance for 
the local-field correction Gk (q) taken from the theory 
of uniform electron gas in the followin g fo rms: 

1) the modified Hubbard correction [20| 


Gk{q) = ^ 


2 g2 -I- fc2 _|_ ^p2 ’ 


( 9 ) 


where ^ is a parameter, the values of which are given 
below; and 

2) the Ichimaru correction 211 
Gkiq) = AQ^ + BQ^ + C + 

2 Q 




4-g" 

X -In 

4g 


2-g 


( 10 ) 
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Fig, 1, Dependences of the dimensionless two-particle 
“density—density” correlation function on the coordinate of elec¬ 
tron 1 normal to the interface, whereas the corresponding co¬ 
ordinate of electron 2 is fixed {z\ = —19.57aB); q = 0.02ag 



Fig, 2. The same as in Fig. 1, but for 2:1 = —2.99aB and q = 
= 0.45a”l 


Table 1. Legends of the figures 


Line type 

Approximation 


Animalu 

— 

Geldart and Vosko = 2) 


Gkiq) = 0 


Hubbard (^ = 1) 

Sham (^ = 1 H- - —) 


Animalu (f = 1 H- - —) 

TrppciB ■' 


where Q = \/+ k'^/pF, and the parameters A, B, 
and C are cumbersome and can be found in work |2l| . 


All numerical calculations were carried out for 
potassium (rs = 4.86 ob). 


Substituting the numerical solution of Eq. ([S]) into 
formula we obtain a two-dimensional Fourier 
transform of the effective potential of electron-elect¬ 
ron interaction. Making the inverse Fourier transfor¬ 
mation with respect to the variable q and taking into 
account that the two-dimensional Fourier transform 
of the effective potential depends only on the absolute 
value of the vector q, we obtain the effective potential 
of electron-electron interaction in the form 

g{r\\,zi,Z2) = ^^e*'i'’ii5(g|zi,Z2) = 

cxD q 

= ^ J(iqqJo{qr\\)g{q\zi,z2), 

0 

where Jo (a:) is the zeroth-order cylindrical Bessel 
function. 

In Figs. 1 to 3, the results of calculations of the 
two-particle “density-density” correlation function of 
electrons, 3Jl(q|2:i, 02 ), obtained by solving the Fred¬ 
holm integral equation of the second kind ([5|) making 
allowance for various local-field corrections and with¬ 
out them (the random-phase approximation) are de¬ 
picted. If one of the electrons is in the metal depth, 
this function is symmetric with respect to the coor¬ 
dinate of a fixed electron directed normally to the in¬ 
terface plane and has a sharp peak, when those coor¬ 
dinates of two electrons coincide (see Fig. I), i.e. the 
electrons correlate with each other and do not feel 
the influence of the surface. If one of the electrons 
approaches the interface, the latter starts to affect 
the two-particle correlation function of electrons. The 
sharp symmetric peak that was observed in Fig. I 
loses its symmetry and broadens: besides electron 
correlations, there emerges an effective repulsion from 
the interface (see Fig. 2). This repulsion results in the 
following. When the electron approaches the interface 
even more, the maximum in the two-particle correla¬ 
tion function of electrons does not occurs, when their 
coordinates coincide, as it was in the metal depth (see 
Fig. 1), but is a little shifted to the left from the in¬ 
terface (in Fig. 3, this maximum is located at about 
-2.7aB). 

In addition. Figs. I to 3 demonstrate that various 
models of the local-field correction do not change the 
behavior of the two-particle correlation function of 
electrons qualitatively, but do it quantitatively. The 
application of the random-phase approximation pro¬ 
duces the smallest deviations in the two-particle cor¬ 
relation function of electrons, whereas the Ichimaru 
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Fig. 3. The same as in Fig. 1, but for 21 = —0.57aB and q = 
= 0.45a-l 



Fig. 4- Dependences of the effective electron-electron inter¬ 
action potential on the distance between the electrons along 
the interface. The electron coordinates 21 and 22 normal to the 
interface are identical and fixed. The case ( 21 , 22 ) — 00 . The 

notation of curves is the same as in Fig. 1 


correction leads to the largest ones. The calculated 
values of two-particle correlation function for elec¬ 
trons with the use of other local-field corrections fall 
within the interval between the values obtained in the 
random-phase approximation and with the use of the 
Ichimaru correction. 

In Fig. 4 to 9, the results of calculations for 
the effective potential of electron-electron interaction 
g{r\\, zi,Z 2 ) obtained for various local-field correction 
models and without them (the random-phase appro¬ 
ximation) are shown. In Fig. 4, the effective poten¬ 
tial of interaction between electrons located in the 
metal depth, i.e. when they do not feel the inter¬ 
face influence, is exhibited. From this figure, one can 
see that making allowance for the local-field correc- 
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Fig. 5. The same as in Fig. 4, but for zi = Z 2 = —Sub 



Fig. 6. The same as in Fig. 4, but for z, = 22 = —Sob 



Fig. 7. The same as in Fig. 4, but for 21 = 22 = —ub 


tion brings about the appearance of a potential well 
at distances from 4.5 to 5 times ob, depending on 
the specific correction model. The depth of this po¬ 
tential well also depends on the local-field correction 
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Fig, 8. Dependences of the effective electron-electron inter¬ 
action potential on the distance from electron 1 to the inter¬ 
face. The normal coordinate of electron 2 is fixed, Z 2 = —5aB, 
ry = = 0. The notation of curves is the same as in Fig. 1 



model. The deepest potential well corresponds to the 
Hubbard correction, and the shallowest one to the 
Sham model. Other examined corrections give rise to 
intermediate depths of the potential well. In addi¬ 
tion, in the random-phase approximation (i.e. when 
the local-field correction is not taken into consider¬ 
ation, the dash-dotted curve), there is no potential 
well at indicated distances in the case of a uniform 
system. At large distances, there emerge Friedel oscil¬ 
lations However, in the case of the effective 

potential calculated in work [1^, they are absent. 


In Fig. 5, the dependence of the effective electron- 
electron interaction potential on the distance between 
the electrons is depicted for the normal coordinates 
of electrons 21=^2 = — Sob- For this distance of elec¬ 
trons from the interface, a considerable deepening of 


the potential well is observed; moreover, it appears 
even in the random-phase approximation. This fact 
originates from a nonmonotonic behavior of the elec¬ 
tron density distribution n(z) near the interface Q; 
namely, here, the plane layers with electron concen¬ 
trations lower and higher than that in the metal depth 
alternate. As a result, the collective effects in the 
electron-enriched layers are more pronounced, and 
the electron screening is stronger. In the electron- 
depleted layers, the situation is opposite: the screen¬ 
ing is weaker, and the repulsion between electrons 
becomes stronger (see Fig. 6). If the electrons come 
nearer to the interface, the repulsion between them 
prevails (Fig. 7). As the coordinates of electrons nor¬ 
mal to the interface grow further, the effective po¬ 
tential of electron-electron interaction tends to the 
Coulomb potential: 


lim g{rii,zi,Z2) = 

Zi,Z2—>CO 


+ (zi - Z2y 


The same behavior is demonstrated in Figs. 8 and 
9. They exhibit the dependence of the effective po¬ 
tential of interaction between the electrons located 
on the same normal to the interface (r|| = 0) on the 
normal coordinate of one of the electrons, regarding 
the other electron to be fixed. The presence of the 
interface results in a nonsymmetric effective poten¬ 
tial of electron-electron interaction with respect to the 
electron coordinate normal to the interface. There are 
more electrons to the left from the fixed one, and the 
screening is stronger; therefore, the potential wells 
and Friedel oscillations are observed; to the right, the 
number of electrons is smaller, so that the potential 
well is either shallower or disappears. 


5. Conclusions 

To summarize, the two-particle correlation function 
of electrons and the effective potential of electron- 
electron interaction have been calculated making 
allowance for various local-held correction models 
known in the theory of uniform electron gas. The 
behavior of the two-particle correlation function of 
electrons depending on the electron coordinates nor¬ 
mal to the interface is studied, as well as the inhu- 
ence of various local-held correction models on it. In 
particular, it is found that the presence of the inter¬ 
face gives rise to an additional effective repulsion of 
the electrons from the interface, the maximum of the 
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two-particle correlation function does not take place 
at the coincidence of electron coordinates, as occurs 
in the metal depth, but is a little shifted from the 
interface. 

The behavior of the effective electron-electron in¬ 
teraction potential depending on the electron coor¬ 
dinates normal to the interface and the distance be¬ 
tween the electrons along the interface is also exam¬ 
ined, as well as the influence of various local-field cor¬ 
rection models on it. The results of our calculations 
demonstrate that, in the near-surface region of the 
metal, there are plane layers, where the behavior of 
the effective potential of electron-electron interaction 
is essentially different. Namely, in some layers, the 
effective potential oscillates and form deep potential 
wells, whereas the wells are shallower or even absent 
in other layers. As a result, the additional mechani¬ 
cal stresses emerge near the metal surface, which can 
provoke the appearance of cracks and other defects. 
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B.M. MapKoeuH, I.M. SadeopunK 

EdiEKTHBHHH HOTEH IT TA.TI 

MDKEJIEKTPOHHOI BBAGMO/fll 

flJIH HAHIBOBMEJKEHOPO EJIEKTPOHHOFO PABY 

3 BPAXYBAHHHM nOHPABKH HA JIOKAJIbHE HOJIE 

P e 3 K) M e 

y po6oTi npoBe^geno uncejibHHH pospaxynoK ecjreKTHBHoro no- 
Tenpiajiy MijKejieKxpoHHOi B3aeMO/i,ii xa .gBOuacxHHKOBOi Kope- 
jiHpiHHo'i 4)yHKU,h ejieKxponiB “rycxHna—rycxnna” fljin naniB- 
obMejKenoro npocxoro Mexajiy 3 ypaxyBannaM nonpaBKH na 
jiOKajiBHe nojie. /tocjii/i,?KeHO BnjmB na hhx njiocKoi' noBepxni 
no/i,ijiy xa pisHHx (jropivi nonpaBKH na jiOKajmHe nojie. HoKaBa- 
HO, ni,o 6ijiH noBepxni no/i,ijiy e objiacxi 3 bijinmoio rjiHbHHOio 
noxeHpiajiBHo'i hmh y ecjreKXHBHOMy noxeHpiajii MijKejieKxpoH- 
HOi B3aeMO/i,ii, nine b rjiHOnni Mexajiy. 
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